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ABSTRACT 

Although all popular approaches to quantum gravity are able to recover the Bekenstein- 
Hawking entropy- area law in the thermodynamic limit, there are significant differences in 
their descriptions of the microstates and in the application of statistics. Therefore they 
can have significantly different phenomenological implications. For example, requiring in- 
distinguishability of the elementary degrees of freedom should lead to changes in the black 
hole's radiative porperties away from the thermodynamic limit and at low temperatures. 
We demonstrate this for the Banados-Teitelboim-Zanelli (BTZ) black hole. The energy 
eigenstates and statistical entropy in the thermodynamic limit of the BTZ black hole were 
obtained earlier by us via symmetry reduced canonical quantum gravity. In that model the 
BTZ black hole behaves as a system of Bosonic mass shells moving in a one dimensional har- 
monic trap. Bose condensation does not occur in the thermodynamic limit but this system 
possesses a finite critical temperature, T c , and exhibits a large condensate fraction below T c 
when the number of shells is finite. 

PACS 04.60.Ds, 04.70.Dy, 03.75.Lm 



1 e-mail address: Cenalo. Vaz@UC.Edu 

2 e-mail address: Rohana.Wijewardhana@UC.Edu 



1 



I. INTRODUCTION 



A solution to the problem of explaining black hole thermodynamics |1|-|3[ from a canonical 
or microcanonical ensemble can be expected to provide important insights into quantum 
gravity, which is considered to be a central problem of theoretical physics. However, there 
is no general consensus on the true nature of the elementary degrees of freedom because 
many apparently different approaches to quantum gravity have all successfully recovered 
the Bekenstein-Hawking entropy. 

The differences between the approaches appear fundamental and range from their intrinsic 
description of the quantum gravity states to the statistics used in counting them. For 
example, in string theory the microstates are dual to weak field D-brane states [§]-|6| , in the 
AdS/CFT approach they are taken to be the states of a particular horizon conformal field 
theory (CFT) [f?HTT| and in loop quantum gravity (LQG) they are represented by punctures 



of a spin network on the event horizon [12-17] . It seems unlikely that these descriptions 
are related in some as yet unknown way because of a more subtle distinction between them, 
which is their use of statistics. In string theory and AdS/CFT the elementary degrees of 
freedom are indistinguishable whereas they must be treated as distinguishable in LQG in 
order to recover the Bekenstein-Hawking law. The statistical properties of the microstates 
can have far ranging phenomenological implications particularly at low temperatures 
and outside the thermodynamic limit. 

In this paper we ask what behavior one may expect outside the thermodynamic limit. 
Such a question can only be answered in a model dependent way, therefore we pose it within 
the context of a model of the quantum black hole that arises from canonical quantum 
gravity applied to spherical gravitational collapse ||19|| . In this model a black hole microstate 
is viewed as a particular distribution of matter shells, which have crossed the horizon, among 
the available black hole energy levels. This is a natural way to realize Bekenstein's original 
ideas in M, from which we quote: "It is then natural to introduce the concept of black hole 
entropy as the measure of the inaccessibility of information (to an exterior observer) as to 
which particular internal configuration of the black hole is actually realized in a given case." 

The simplest example is provided by the static Banados-Teitelboim-Zanelli (BTZ) black 
hole p0[ . A symmetry reduced canonical quantization of this black hole yields Hawking 
radiation |2l| and an equispaced black hole mass spectrum. Furthermore, counting the black 
hole' s microstates in the canonical ensemble reproduces the Bekenstein-Hawking entropy in 
the thermodynamic limit, provided that the counting proceeds using Bose statistics p2[ . The 
thermodynamic limit is the limit in which both the system size and energy are simultaneously 
taken to be very large. For the BTZ black hole this means that the temperature is high. Our 
interest here is to examine the same ensemble at low temperatures and for a finite number 
of matter shells. We will show that Bose-Einstein condensation occurs and the BTZ black 
hole turns "cold" once condensation sets in. "Cold" black holes are of course stable against 
Hawking evaporation. 

One may criticize our choice of the 2+1 dimensional black hole, even though it has the 
virtue of being simple, because it is topological and therefore unique in many respects. 
For instance when canonical quantum gravity is applied to the Schwarzschild or higher 
dimensional black holes in the absence of a cosmological constant the microstates are counted 
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in an "area ensemble" using Boltzmann statistics. However the situation is different for 
"large" Anti-de Sitter (AdS) black holes, which are black holes whose horizon radius is 
much larger than the AdS length. In a particular limit, their mass spectra are identical 
to the spectrum of the BTZ black hole modulo dimension dependent constant factors and 



Bose statistics must be applied to recover the Bekenstein-Hawking area law p3fl . Moreover, 



many higher dimensional black holes from string theory have the form BTZ x Ai, where Ai 
is a simple manifold. Their thermodynamic properties can be recovered directly from the 
thermodynamic properties of the BTZ black hole p4| |25|| . Bose condensation is a possibility 
whenever indistinguishability is an essential feature of the counting of microstates. It may 
thus occur for realistic black holes in higher dimensions, in which case one can expect stable 
black hole remnants to contribute significantly to the Dark Matter content of the universe. 

This paper is organized as follows. In section II we review the canonical quantization of 
the BTZ black hole and its statistical thermodynamics in the canonical ensemble, as given 



in We then go outside the thermodynamic limit in section III, treating the black hole 
as a system of finite size. With the black hole spectrum, which we derive in section II, we 
show that Bose condensation can occur. We go on to determine the critical temperature and 
the condensate fraction as well as the entropy and argue that the area law strictly holds only 
in the thermodynamic limit but breaks down outside it and below the critical temperature. 
We conclude in section IV with some comments. 



II. STATISTICAL THERMODYNAMICS OF THE BTZ BLACK HOLE 

The static BTZ black hole may be viewed as the end state of the collapse of inhomoge- 
neous dust |26|, p7| , which is described by the LeMaitre-Tolman-Bondi family of solutions 
of Einstein's equations with a negative cosmological constant. The solutions are character- 
ized by two arbitrary functions, viz., the "mass function", F(p) = 4GM(p), representing 
the initial mass distribution of the dust ball and the energy function, E(p) representing its 
initial energy distribution. The classical solutions are given by the circularly symmetric line 
element 

ds 2 = dr 2 - ^r^dp 2 ~ RV, (1) 

where r is the dust proper time and p labels dust shells or curvature radius R(t, p). Einstein's 
equations can be integrated to give the energy density of the dust 

and a dynamical equation for the evolution of the area radius, 

(d T R) 2 = 2E- AR 2 (3) 

in terms of two arbitrary integration functions of the shell index coordinate, p. The last 
determines the area radius of shells to be 



R ( T 'P) = \/¥ sin ( -v/Ar + sin" 1 jA p ) } ( 4 ) 
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where we have used a freedom in the scaling of the shell index p by setting R(0, p) = p and 
the fact that collapse solutions satisfy d T R < 0. The integration function F(p) represents 
the initial mass distribution of the collapsing dust ball and the function E(p) represents the 
initial energy distribution; one finds 

F(p)= [ P e(0,p)pdp 
Jo 

E(r) = 1 -[(d T R)l =0 + AR 2 ]. (5) 

directly from @ and (||) respectively. 

The general circularly symmetric Arnowitt-Deser-Misner (ADM) line element, 

ds 2 = N 2 dt 2 - L 2 (dr - N r dt) 2 - R 2 d<p 2 (6) 

can be embedded into the metric in ([!]). After a series of transformations described in detail 



in 



2lH this leads to a canonical description of the classical black hole in terms of the dust 



proper time, r(r), the area radius, R(r), the mass density defined via 



M, 



Hr) = I dr> T(r), (7) 

* Jo 

where M is an arbitrary constant contributed by the boundary at the origin, and their 
conjugate momenta, P T (r), Pr{t) and Prir) respectively. The effective constraints of the 
gravity-dust system are then obtained in the form 

H r = t'P t + R'P R - TP^tt 

n = P 2 + TP 2 -y^0, (8) 

where J 7 = AR 2 — F and the prime refers to derivatives with respect to the ADM label 
coordninate r. 

Applying Dirac's formal quantization to the above system of constraints, we replace 
the momenta by functional derivatives with respect to their corresponding configuration 
variables. To encapsulate the factor ordering ambiguities at the formal level of the Wheeler- 
DeWitt equation we introduce factors of (5(0) into the resulting functional Schroedinger 
equation, writing the quantum Hamiltonian constraint as 



r 



wh[t, r, r] 



r2 r2 r p2' 



*[r,fl,r]=0, (9) 



where A(R,F) and B(R,F) are smooth functions. The continuum limit of the wave- 
functional is taken to be of the form 



V[t,R,T] = exp 



i / drT(r)W{r(r),R(r),F{r)) 



(10) 



which formally obeys the momentum constraint provided that W(r, R, F) has no explicit 
dependence on the label coordinate r. We then put @ and (|10D on a lattice and take the 



4 



continuum limit p8| . This corresponds to the choice of regularization. Consistency of the 
lattice regularization fixes the factor ordering via a set of three equations |29j for each lattice 
site, j, 



dR 



1 



o. 



d- 2 



j 



dR 2 



+ A 3 



OR, 



B 4 = 0. 



(11) 



Thus B(R,F) is constrained to be identically vanishing and A(R,F) and the lattice wave- 
functions are obtained by solving the first two equations. Furthermore, hermiticity of the 
Hamiltonian constraint requires that the Hilbert space measure m(R, F) is determined from 
A(R, F) according to 

A 3 = |^|9 R Tn(m,|^-|). (12) 

Thus, the system is completely solved once regularized on the lattice. 

Unfortunately, the situation becomes more complicated when the mass density function 
is distributional, as it is in the description of the final state black hole, for then the system 
naturally collapses into a countable product of wave functions, the functional differential 
equations become ordinary partial differential equations and no regularization is required. 
It follows that no further conditions are available, leaving the factor ordering ambiguity in 
the form of the unknown functions A and B unresolved. 

This is exemplified by the BTZ black hole whose metric is of the form 



ds 2 = - (AR 2 - 8GM) dT 2 + 



dR 2 



(AR 2 - 8GM) 



+ R 2 dip 2 



(13) 



where M is the black hole mass parameter. This is a special solution of (P in which the 
mass function is taken to be constant, F = 4GM, for p > and the energy function is given 
by 2E = 1 + 8GM, again for p > 0. The metric in (0) can be brought to the static form in 
([13|) by the transformations R = R(t, p) as given in @ and 



T = r+ dR 



.V l + 8GM - AR 2 
AR 2 - 8GM 

for the Killing time, T. Within the framework of the canonical theory we take 

F(r) = AGM + AGeQ(r) = AGM, 



(14) 



(15) 



where e represents the mass of a shell at r = and where G is the Heaviside function. 
Likewise, the energy function should be given by 

E{r) = l - [1 + 8GM + 8GeO(r)] . (16) 

We see that the mass function in fli"5|) yields a mass density that is the 5— distribution 

r(r) =AGe5{r), (17) 



5 



and the wave-functional in (|I]) turns into the wave-function, 



where r = r(0), R 



R(0) and F 
d 2 



J™ drT(r)W(T(r),R(r),F(r)) 



^ieW{T,R,F) 



F(0). The Wheeler-DeWitt equation becomes 



Or 2 



+J: m +A m +B 



ieW(r,R,F) 







(18) 



(19) 



and we note that F) and B(R, F) remain undetermined, as does the measure m(R, F), 
with B(R,F) playing the role of an external potential. 

Fortunately there is a way out of this difficulty which is provided by a detailed examina- 
tion of the Hawking radiation from the black hole within the context of this midisuperspace 
model. Because the Wheeler DeWitt equation is second order in the dust proper time we 
can define both positive and negative frequency states with respect to r. Therefore, taking 
the collapsing dust as a small perturbation around a pre-existing, massive black hole, so 
that the role of the quantum matter in Hawking's original background field derivation of 
black hole radiance |J is played by the dust, the system (|ll|) can be shown to yield Hawking 
radiation at the Hawking temperature [21]. However, an important point is that obtaining 



the correct Planckian distribution in the near horizon limit requires a particular choice of 
measure, one appropriate to the massive black hole because the dust itself is considered a 
perturbation on the black hole background. The measure turns out to be the one obtained 
from the deWitt supermetric, which, from (H), is 



Thus 



lab 



m 



1 

0^ 



(20) 



(21) 



and putting this together with the hermiticity requirement in (12) determines A(R, F). Now 
since the black hole ends up being a single shell in this simple quantum mechanical model, 
it is reasonable to take the external potential, B, to be vanishing. Thus ( |I9"D reduces to 
the free Klein-Gordon equation which is hyperbolic in the interior and elliptic elsewhere. A 
stationary state solution this equation was shown in |22| to yield the mass levels 



3 + 



(22) 



where I 2 = — A -1 is the AdS length. 

The macroscopic black hole is defined as the end state of the collapse of many, say 
M = ^2jAfj, matter shells, Mj of which occupy level j of (22). Thus, when a boundary 
contribution from the origin is included [Ef| |27j 3 , the total black hole mass is given as 



M = Mi 



h 



+£7 U + 



(23) 



3 In more than three dimensions this contribution from the origin is usually set to zero, otherwise it would 
represent a singular initial configuration. In three dimensions a non- vanishing Mq is necessary to allow for 
a velocity profile that vanishes at the origin. This does not lead to singular initial data and the presence 
of Mq does not lead to a singular initial configuration. 
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From this point of view, the elementary degrees of freedom are bosonic mass shells and a 
black hole microstate is a particular distribution of M shells between the levels in (p2|) . 

The statistical thermodynamics, specifically the H — > oo limit, of this model can be 
treated in the canonical ensemble and is captured by the partition function 



i.v 



exp 



-(3 (m +j2^^J 



(24) 



where A/} represents the number of shells excited to level j, with mass Ej, and 
g(Af\, . . . ,Afj, . . .) is the degeneracy of states, which we take to be unity so as to imple- 
ment Bose statistics. The canonical entropy is then obtained from 



S am =\pM + ]nZ(P)] lli 



M=-dlnZ/d/3 ' 



(25) 



where M is the average energy in the canonical ensemble, which we associate with the black 
hole mass. When the system size is taken to be infinite, the partition function 



3=0 



(26) 



is easily evaluated. It may be rewritten, by exploiting the well known duality in as 

1 /tt 2 ! , Ph\( SlM n 1 



Z((3h/2l) = —e ^ ~«>[Z{A<K 2 l/ph)] 

\/2 



-i 



(27) 



This links the high temperature behavior of our system to its low temperature dynamics. If 
one assigns a value Aq to the ground state energy of the system, taking 



Z(A7T 2 l//3h) w e 
then it is then a straightforward exercise to show that 



4n\ a 



■eft' 



IM 

6h 



where 



Ceff 



l-^(Mo-A Q ) 
n 



(28) 



(29) 



(30) 



depends on two undetermined parameters, the vacuum energy A and the boundary contri- 
bution Mq. The vacuum energy depends on what one takes to be the ground state of the 
BTZ solution, for example the choice A = —1/SG corresponds to the choice of pure AdS 3 , 



dR z + R 2 d<p 2 



(31) 



for the ground state |JT]. On the other hand, the presence of the arbitrary contribution, 
Mq, from the boundary at the origin must be fixed by some other means, for example by 
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comparing ( p9|) to the Bekenstein-Hawking entropy of the BTZ black hole. The Bekenstein- 
Hawking entropy is just 

A h nR h Til [2M 
SB - u -4G-2Gh-J\l-G> (32) 
where Ah is the horizon area and Rh its radius. This in turn implies that 

Ccff= 2^' (33) 

and therefore fixes M — A . 



III. THE GRAND CANONICAL ENSEMBLE 

We now consider a black hole made of a finite number M = . A/} shells. It is standard 
practice to then begin with the grand partition function, 

S(e, f3) = e~ m ° Y[(l - e ~^'^)-\ (34) 

3 

where we introduced a chemical potential /i, which is determined from the constraint 

£<A/->=A/- (35) 

3 

and where Ej is the energy of the level j. Introducing A = e^ M ~ e °) < 1 (A -1 is the fugacity) 
and Aej = Ej — Eq we then have 

^fl-^IlOT (») 

3 

The average level occupation number is computed directly from E((3,fi), 

W> = -y^ - (37) 



so our condition for determining the chemical potential is 

3=0 3=1 



00 1 A °° A 

* = E ^ _ A = + E ^ _ A • (38) 



The first term on the right represents the ground state occupancy (henceforth, average values 
will be understood and the angular brackets will be omitted), 

No = y^a' ( 39 ) 
and the second represents the number of excited shells, 

N - = J2^—-y (40) 
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The condensate fraction is defined as the ratio of the ground state occupancy to the total 
number of shells, fc = N /Af. 

We will first determine fc as a function of the temperature. For the spectrum of the BTZ 
black hole in ([22]), Aej = hj/l and therefore, replacing the last sum in (p8|) by its power 
series expansion, we can express the total number of shells as 

M = + A r e-^/*. (41) 

j=l r=l 

The most direct way to evaluate the double sum in fl4l|) is to employ the Mellin-Barnes 
representation of the exponential function, 



i pr~{-ioo 

— : / dt T(t)a- 1 (42) 

271"^ Jr—ioo 



where r G M and Re(a) > 0, to re express it as an integral over the complex plane. This 
gives 

A 1 r +<~ r(f)C(«)Li,(A) 



M —x + ^u dt ^mt^- <43) 

where ((t) is the zeta function and Lii(A) is the polylogarithm. For A < 1 the integral is 
governed by the simple pole of the (— function at t = 1, therefore 

Equation ( PD can be solved for A in terms of the Lambert function, W(x), 

A <^> = 1 " wmielUm) («) 

and determines the ground state occupancy as well as the condensate fraction, 

f c %N)-jf- - mK[l ^ (46) 

in terms of the temperature. It is more illuminating, however, to express the condensate frac- 
tion as a function of the dimensionless ratio, x = T/T c , where T c is the critical temperature 
at which all the shells are excited. 

To calculate the critical temperature, T c , it is necessary to set A ~ 1 and iV cx = Af, the 
total number of shells. Once again employing the Mellin-Barnes transformation, this time 



replacing the polylogarithm function Li f (A) |a^i by the (— function, we find [[32 



N i p mew 

but the integral now admits a double pole at t = 1 via the (— function and therefore 

( r(t) y / 



=-s in(ew) - (48) 
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where 7 is the Euler-Mascheroni constant. Solving for the inverse critical temperature, we 
determine 

Wl = M , (49) 



which is the exact version of the result in The Lambert function increases linearly for 



small values of its argument and slower than the log function for larger values. Thus, for a 
few shells the critical temperature is more or less constant whereas it increases roughly as 
M I In M when M grows large. This is a well known property of one dimensional harmonic 
traps. When the black hole is composed of many shells, most of them will in fact lie in the 
ground state. However, for a small number of shells the ground state occupancy can still be 
a large percentage of the total number of shells in a significant interval of temperatures. 



Using (|49 ) to re-express the condensate fraction in terms of the variable x = T/T c , we 



easily arrive at 

n / »rs x (■V+i)W JV 1 

k(x ' N) = m w{ a^ cT) - JT (50) 

where we have defined WV = W(e~ 7 A/'). At high temperatures, 



whereas, at low temperatures, 



, x-+o 1 InjV ... 2 . 

/c ~ 1 - Tjn-z + °( x ) (52) 



iV 



Figure [T] shows the behavior of fc as a function of x for A/" ranging from one hundred to 
one hundred thousand shells. It is small, increasing slowly to the left when x 3> 1 but 
significantly faster when 1 < 1 and approaching unity in the limit 0. Again, 

the slope of the condensate fraction, dfc/dx, approaches In AT / Wn for small x and it is 
vanishingly small for large x. Figure |2] shows this slope for the same range of M. Both the 
analytical approximations above as well as the figures indicate that increasing the number of 
shells saturates the slope of the condensate fraction. This behavior suggests that the critical 
temperature may be thought of as the condensation temperature. 

The average energy of the system, which is associated with the black hole mass, can be 
expressed in terms of the occupancy of the levels in the usual way and the sum rewritten in 
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TIT, 



FIG. 3: C/N as a function of T/T c 
the integral representation we have used so far; one gets 

h/l 



M{P,Af) = M +ATe + 



^ r(t)Li t (A)C(t-i; 

2™ ./ T _ ioo (/3/i/0* 



(53) 



For A < 1 there is only a simple pole coming from the (^—function at t = 2, so the black hole 
mass can be given in terms of the temperature by 



M(/3,AT) = M +Afe + 



Li 2 (A) 

(3 2 h/i 



and the heat capacity is 



(54) 



(55) 



AT 



where the prime on the polylogarithm represents a derivative with respect to A. Its behavior 
is best understood in terms of the dimensionless parameter x = T/T c , 



Wat \W n \oxJ n 



(56) 



The specific heat, C/N", is displayed as a function of x = T/T c in figure ^] for ranging 
from one hundred to one hundred thousand shells. There is no true phase transition in this 
system. At low temperatures, the fugacity in fl45 ) behaves roughly as 1 — A/" -1 and the 
second term in the expression for the heat capacity in flSEQ becomes vanishingly small. Thus 
in the limit as (3 — > oo, using Li 2 (l) = 7r 2 /6 we find 



C 



7l 2 kT 



N Vsfh/l' 

(see |M| ^]). On the other hand, at high temperatures the fugacity behaves as 

A = 1 - e-W 

and the specific heat approaches unity in the limit as — > 0. 
We now turn to the entropy, which may be computed from 



(57) 



(58) 



S = \nE + 0U - fiPJ\f. 



(59) 
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FIG. 4: S/N as a function of T/T c . 



Expanding In H in a power series and applying the transformation in (03) we arrive at 



lnS = -/3M -ln(l-A) + ^^ (60) 



and combine this expression with (p4|). Then the entropy function can be given in terms of 
the ratio x = T/T c as 

S(P,Af) = ^Li 2 (A) - ln(l - A) - ATln A. (61) 

VVN 

The specific entropy is shown in figure |] for TV" ranging from one hundred to one hundred 
thousand shells. It is seen to decrease with increasing N '. In the limit of low temperatures, 
i.e., as x — > 0, the leading behavior of the entropy is as 

This is not vanishing, but the specific entropy, S/Af, does vanish in the thermodynamic 
limit as the temperature approaches zero and so it is consistent with a generalized version 
of the third law. If we now consider this function in terms of its natural variables, M and 
M we find from (54) that in this same limit 



vr ^ /6(M-M -A/£o) ; (g3) 



fih/i y h/i 

showing that the black hole temperature vanishes as the black hole mass approaches Mq + 
Neq, at which point the entropy approaches lnA^. In fact, written in terms of the black hole 
mass and the number of shells, the entropy is 



s ~W 6(M ~yr^ o)+inAf+i - (64) 

Small x is the regime in which Bose condensation plays a dominant role. At temperatures 
well below the critical temperature, where a large fraction of the shells are in the ground 
state, the first term on the right becomes negligible and the entropy is dominated by its 
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AjV„./JV 




FIG. 5: AN ex /Af as a function of T/T c . 



logarithmic behavior. Thus even though a quantization of Einstein's gravity yields the usual 
black hole results in the thermodynamic limit and at high temperatures, its behavior outside 
this regime may present interesting features peculiar to quantum statistics. In particular, 
we see that BTZ black holes turn cold once condensation sets in. This is the "small" black 
hole or short distance regime in 2+1 dimensions. 

Finally, for completeness, we address the mean square fluctuations. We are interested in 
evaluating 

AA, « = -g E It = E £ rX ' e ~^"- (65) 

^ j=l i j=l r=l 

The sums can be evaluated with the help of the representation in ( pE2]) and one obtains 

2 = Li 2 (A) 

ex (/3fr//) 2 1 ; 

The fluctuations are plotted as a function of x — T/T c in figure |^ for one hundred, one 
thousand and one hundred thousand shells. The uppermost curve represents M ranging 
from one hundred to one hundred thousand shells, and the fluctuations are seen to decrease 
with N '. Near the critical temperature, they fall off as 1/ In M . 



IV. DISCUSSION 

Although well adapted to problems possessing high symmetry such as in cosmology and 
spherical collapse, the canonical approach suffers from several ambiguities. We have dealt 
with them explicitly before arriving at the results quoted in the introduction, but only within 
the context of the specific class of models with which we work. Thus canonical quantization 
cannot be expected to yield the final theory of quantum gravity as it stands, but its results 
should be taken seriously for models in which the ambiguities can be resolved. One advantage 
of our approach to the quantum properties of black holes is that the fundamental degrees 
of freedom and the black hole microstates are given a transparent meaning. Another is that 
the same picture holds for higher dimensional black holes in Einstein's general relativity 
both with and without a cosmological constant. 

It seems remarkable at first sight that the spherical midisuperspace models of collapse 
should yield the Bekenstein-Hawking entropy of the black hole, considering that all non- 
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spherically symmetric states have been eliminated by the symmetry reduction and therefore 
do not contribute to the entropy. We cannot definitively tell why this is so, but we speculate 
that non-spherically symmetric states may not describe a regular horizon (with constant 
surface gravity) and so do not contribute to the entropy of a static black hole. However, 
this conjecture requires further clarification and the question seems worthy of further inves- 
tigation. 

The point we have made in this paper is that the statistics applied in counting black hole 
microstates, together with the specific models of the microstates, can play a crucial role in 
the behavior of a black hole at low temperatures and outside the thermodynamic limit. We 
have demonstrated this explicitly in the case of the BTZ black hole for which Bose statistics 
must be applied to correctly recover its thermodynamic properties. We have argued that 
the BTZ black hole condenses in the sense that there is a critical temperature below which 
a significant fraction of the black hole mass is in the ground state. 

Our results can be summarized as follows: (i) the condensation temperature grows almost 
linearly with the number of shells and can be quite large, (ii) the specific heat grows linearly 
below the condensation temperature but slowly above it, approaching unity in the limit of 
very high temperatures and (iii) the Bekenstein-Hawking area law is obeyed only in the 
thermodynamic limit, outside of which the entropy deviates from the area law, approaching 
In A/* in the limit as T — > 0, which occurs at finite mass. According to the results of p6| |, 
this would signal a breakdown of general relativity at low enough temperatures. This seems 
to be an interesting direction for future research. 

Of considerable phenomenological importance of course is the question of whether these 
results are reproduced by realistic black holes in higher dimensions and precisely how. If 
they are, then Bose condensed primordial black holes would make natural Dark Matter 
candidates. We will examine the consequences of the Bose condensation described in this 
paper for gravity at short distances in a future publication. 
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